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Proximate transition temperatures amplify linear magnetoelectric coupling in 
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We report a giant linear magnetoelectric coupling in strained BiMnOs thin films in which the 
disorder associated with an islanded morphology gives rise to extrinsic relaxor ferroelectricity that 
is not present in bulk centrosymmetric ferromagnetic crystalline BiMnOs. Strain associated with 
the disorder is treated as a local variable which couples to the two ferroic order parameters, magne¬ 
tization M and polarization P. A straightforward “gas under a piston” thermodynamic treatment 
explains the observed correlated temperature dependencies of the product of susceptibilities and 
the magnetoelectric coefficient together with the enhancement of the coupling by the proximity 
of the ferroic transition temperatures close to the relaxor freezing temperature. Our interpreta¬ 
tion is based on a trilinear coupling term in the free energy of the form L ■ (P x M) where L is 
a hidden antiferromagnetic order parameter, previously postulated by theory for BiMnOs. This 
phenomenological invariant not only preserves inversion and time reversal symmetry of the strain- 
induced interactions but also explains the pronounced linear magnetoelectric coupling without using 
the more conventional higher order biquadratic interaction proportional to (P ■ M) 2 . 

PACS numbers: 75.85.+t 77.55.Nv 77.80.Jk 


I. INTRODUCTION 

Magnetoelectric (ME) coupling, the induction of elec¬ 
tric (magnetic) polarization via magnetic (electric) fields, 
represents an area of intense research due to both its 
complex physical origins and its potential technological 
applications M- According to general thermodynamic 
considerations, multiferroics - materials which simultane¬ 
ously possess spontaneous magnetic and electric polariza¬ 
tions - have long been predicted as prime candidates for 
enhanced ME coupling [3H f|. Additionally, strong limi¬ 
tations on ME coupling, such as mutually exclusive re¬ 
quirements in electronic configurations for each ferroic 
order, have been identifiedQ. Accordingly, understand¬ 
ing global energetic constraints on ME coupling - and 
ways to relax these constraints in multiferroic systems - 
is of fundamental importance for this class of materials to 
reach their potential. In particular, the use of concomi¬ 
tant ferroic transition temperatures to amplify ME cou¬ 
pling has been proposed(6j,[8], but an experimental study 
providing insight into the thermodynamic constraints on 
such an enhancement has yet to be presented. 

Crystalline BiMnOs (BMO) is recognized as the only 
insulating perovskite that is strongly ferromagnetic. 
Early indications that BMO should also be a displacive 
ferroelectric earned it the distinction of being called 
the “hydrogen atom” of multiferroics Orbital or¬ 
dering of Mn spins was determined to be responsible 
for ferromagnetism p~ol| . and a structural distortion in¬ 
duced by the lone pair of the Bi ions was believed 
to produce spontaneous electric polarizationfuQ]. How- 
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ever, subsequent neutron scattering and density func¬ 
tional theory (DFT) investigations have contested the 
non-centrosymmetric classification of BiMnOs, casting 
doubt over its intrinsic ferroelectric properties 11 2| 1 1 -il | . 

Despite these claims of a centrosymmetric crystal 
structure for bulk BMO, reports of ferroelectricity in thin 
films of BiMnOs remain [15il21j . suggesting alternate ex¬ 
trinsic mechanisms are present. In a comprehensive re¬ 
view of polar and nonpolar phases of BiMnOs Belikf22| 
has pointed out that any particular BMO thin-film sam¬ 
ple can have strain-sensitive structural and compositional 
modifications that are strongly correlated with the ap¬ 
pearance of ferroelectricity. Importantly, the film resis¬ 
tivity must be sufficiently high to insure that leakage cur¬ 
rents do not reduce the charge separation associated with 
remnant electric dipole moments. 

We show that in our pulsed laser deposited multifer¬ 
roic BiMnOs thin films the operating extrinsic mecha¬ 
nism that gives rise to a pronounced and surprisingly 
large linear magnetoelectric effect is strain disorder. We 
observe temperature-dependent transitions of the mag¬ 
netic (M) and ferroelectric (P) order parameters which 
are in close proximity to each other and give rise to a 
strongly enhanced magnetoelectric coupling in the vicin¬ 
ity of the transition temperatures. The interacting ferroic 
order parameters and their contribution to a large mag¬ 
netoelectric effect are described using a strain-dependent 
trilinear interaction term in the free energy of the form 
L ■ (P x M). To preserve inversion and time reversal 
symmetry, a staggered order parameter L representing 
an antiferromagnetic background is used. The trilin¬ 
ear form of the interaction originally discussed by Fen- 
nie [23| is motivated by a Dzyaloshinsky-Moriya[24l [25| 
type antisymmetric interaction. Recent theoretical treat¬ 
ments suggest that “hidden” antiferromagnetism must be 
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FIG. 1: Temperature dependence of remanant polarization 
Pr (solid black) and magnetization M (dashed red) demon¬ 
strating the near overlap of the ferroelectric and magnetic 
transitions. The vertical arrows show the temperatures where 
the maxima in the slopes occur. 

present in BiMn0 3 [Il,[27|. Using the more conventional 
higher-order biquadratic interaction, (P ■ M) 2 , which is 
also symmetry allowed, does not explain the data as will 
be explained below. 


II. METHODS 

Our multiferroic samples are 60 nm-thick thin films of 
stoichiometric BiMn0 3 grown using pulsed laser deposi¬ 
tion on (001) SrTi0 3 (STO) substrates. X-ray charac¬ 
terization showed str ong peaks demonstrating coherent 
local order (see Ref. |2yl for stoichiometric analysis and 
growth parameters). The films grow with a 001 orien¬ 
tation (pseudo-cubic notation) on STO substrates with 
a compressive strain due to lattice mismatch of 0.77%. 
Optimized tuning of deposition rate, substrate temper¬ 
ature, target stoichiometry, oxygen partial pressure and 
post deposition cooling rates leads to films with no dis¬ 
cernible impurity phases and resistivities higher than 10 6 
f l cm in the temperature range where all measurements 
were taken [20| . 

The surface topography of the films imaged using 
tapping-mode atomic force microscopy (AFM) shows 
a three-dimensional island growth mode with an rms 
roughness of approximately 10 nrn. Accordingly, nonuni¬ 
form strain with high values of strain at the island edges 
is expected [2^, [29|. This scenario is consistent with 
temperature-dependent magnetization measurements ac¬ 
quired with a Quantum Design SQUID magnetometer 
using an in-plane field of 500 Oc|20|. The magnetiza¬ 
tion shown in Fig. [I] (dashed red line) shows an onset 
near 85 K and a saturated moment of approximately 



FIG. 2: Polarization loops at 90 K acquired on the time scales 
indicated in the legend show relaxor ferroelectric behavior. 
On the longer time scales the remnant polarization has more 
time to interact with the thermal background and decrease in 
magnitude. 


1 /is/Mn. The transition and saturated moment in our 
thin films is lower than the transition ( T c « 105 K) 
and saturated moment (3.6 /ig/Mn) obtained in poly- 
crystalline samples [Hi [3til |. Although the presence of 
Bi vacancies may contribute to reduced magnetism in 
BMOdl, we attribute our reductions of transition tem¬ 
perature and saturated moment primarily to the sub¬ 
strate induced strain and the strain disorder associated 
with the islanded morphology. 

The remanent ferroelectric polarization was measured 
using a Radiant Technologies ferroelectric tester and a 
previously-described interdigital capacitor geometry [20|. 
The polarization in the remanent hysteresis loop is cal¬ 
culated by isolating the transferred charge due solely to 
domain-switching. 

The ferroelectric polarization loops shown in Fig. [2] 
shrink in area as the time (indicated in the legend) to 
complete a loop increases. This behavior is a signature 
of relaxor ferroelectricity in which the charge separation 
associated with induced polarizations decays with time 
because of interactions with the phonons of the ther¬ 
mal background. The remanant polarization Pr for each 
loop is taken as the polarization at an applied electric 
field E = 0 where the horizontal portion of each loop 
crosses the vertical axis. Dielectric characterizations of 
these same films [fill] show two thermally-activated loss 
peaks that are signatures of relaxor ferroelectricity with 
the real part of the dielectric constant becoming fre¬ 
quency inde pen dent and diverging at a relaxor freezing 
temperatureI32i| Tp near 70 K. For T > Tp any rema- 
nence in the sample will decay to zero on laboratory time 
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scales. 

Since Bi-doped SrTiC >3 has been reported to be a re- 
laxor ferroelectric (3l|. we necessarily should be concerned 
with possible ferroelectric contributions from Bi dop¬ 
ing of the SrTiC >3 substrate during the high tempera¬ 
ture growth (632°C) [20j of the BiMnC >3 films. This 
source of ferroelectricity was ruled out via two methods. 
First a BiMnC >3 film was grown and then etched with a 
KI solution. After depositing interdigital electrodes di¬ 
rectly on the exposed substrate, remanent polarization 
measurements revealed distorted hysteresis loops which 
were more than 100 x smaller than those measured for 
BiMnC> 3 . In addition, dielectric measurements did not 
display the relaxations observed in the BiMnOs films. 
In our second method, BiMnOs films were grown on al¬ 
ternate substrates (NdGaOs and SrLaGa 04 ) which have 
different composition from that of SrTiC> 3 . Our measure¬ 
ments showed reduced but open remanent polarization 
loops as expected since the mismatch strain is less. Ac¬ 
cordingly, Bi-doped SrTiC >3 was definitely ruled out as 
the source of relaxor ferroelectricity reported here. 


III. RESULTS AND DISCUSSION 

The proximity of the magnetic and ferroelectric transi¬ 
tions as measured by the overlap of the rapid changes in 
the temperature-dependent magnetization M{T ) (the or¬ 
der parameter for the magnetic transition) and the rem- 
anant polarization Pr(T) is shown in Fig.[Q The relaxor 
(diffuse) nature of the ferroelectric transition shown in 
Fig. [5] implies the presence of disorder-induced precursor 
polar microregions which with decreasing temperature 
grow in size and fluctuate on longer time scales [34|. In 
the temperature regions of Fig. [T| where there is a pro¬ 
nounced increase in the magnitude of Pr with decreasing 
temperature, the polarizations become increasingly sta¬ 
ble in time. The question of whether the proximity of 
the two transitions is accidental will be addressed below. 

As shown in Fig. the polarization loops measured 
at the temperatures indicated in the legend steadily 
increase in area as the temperature is reduced. The 
remanent-polarization ( Pr ) evaluated at the electric field 
E = 0 appears near 110 K and increases slowly to a 
large value of 23/zC/cm 2 at 5K. Application of an ex¬ 
ternal tensile strain on the order of 0.01% using a three- 
point “beam-bending” technique increases Pr at 5 K by 
50% (31.]. Our experimental observations of relaxor fer¬ 
roelectricity, which requires the presence of disorder- 
induced polar micro regions [34j. together with the above 
mentioned sensitivity of Pr to externally applied strain 
strongly suggests that disordered strain regions, which 
are not present in crystalline bulk BMO, are responsible 
for the robust ferroelectricity observed in thin-film BMO. 

In addition to the robust extrinsic ferroelectricity, we 
also observe a strong magnetoelectric (ME) coupling 
which is especially pronounced in the region where the 
transitions in Pr and M overlap. The lower panel of 
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FIG. 3: Temperature and magnetic field dependence of elec¬ 
tric polarization loops. (a) Hysteretic polarization loops 
at the indicated temperatures span a temperature range of 
« 100 K with a maximum remanent polarization, Pr, of 
«23 fi C/cm 2 at 5K. (b) A magnetic field of 7T (red) is 
shown to decrease the FE polarization by « 10% (red) from 
its value in zero field (black). 


Fig. [3] shows how the application of a 7 T magnetic 
field shrinks the area of the polarization loop and re¬ 
duces Pr by S3l0%. The reduction of Pr is linearly 
dependent on the applied magnetic field B = fioH as 
shown in the selected curves of Fig. [4] for the temper¬ 
atures indicated in the legend. We find that the ME 
effect, characterized by a coupling coefficient defined by 
the slope a = —dPn/d^oH), is negative and linear at 
all temperatures. The coupling coefficient is found to be 
independent of field direction and is quite large, reach¬ 
ing a value of -0.1 /iC/cm 2 T (-1.25 ns/m in SI units) at 
T = 65 K. This value is approximately 35 times larger 
than the current record for single-phase linear ME cou¬ 
pling of 36.7 ps/m found in TbPC^P, HI]. Coupling 
larger than that presented here has been reported, but 
strong non-linearities have made such couplings difficult 
to interpret thermodynamically [H, 37]. In addition, a 
averages out the effect of disorder in our sample, since 
the measurement length scales for both M and P are 
much larger (an order of magnitude or more) than the 
typical island size and separation [2C|. The relevance of 
this point will be discussed in section C. 
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FIG. 4: The remnant polarization Pa decreases at all tem¬ 
peratures with a linear dependence on magnetic field. The 
magnetoelectric coupling constants a(T) for the three selected 
temperatures in the legend are extracted from the slopes of 
each of these curves. The maximum d(T)with a value of -0.1 
pC/cm 2 T (-1.25 ns/m) occurs near 55 K. 


A. Gas under a piston thermodynamics 

We have chosen the relevant thermodynamic variables 
describing our multiferroic system to be the order pa¬ 
rameters M and P rather than the magnetic induction 
B and the electric displacement D, since the “medium 
in an external field” is analogous to the “gas under a pis¬ 
ton” system^. By making the analogy to a perfect gas 
described by the state variables p (pressure), V (volume) 
and T (temperature), we incorporate the intensive vari¬ 
ables, M and P, into the free energy bv replacing the 
work term pdV with the substitutions [381. |39| p —>• H, 
V — p>qMV = —p>oXmHV for the magnetic energy 

and p —1 E, V —> —PV = —eoXpEV for the electric 
energy. SI units are used throughout. 

These substitutions become particularly useful in find¬ 
ing expressions for the electric and magnetic susceptibil¬ 
ities beginning with the well-known expression for the 
heat capacity difference C p — Cy of a pVT system , 

C p -C v = ~T{dV/dT\ p f/(dV/dp) T , (1) 

which reduces to C p — Cy = R for one mole of a per¬ 
fect gas with R equal to the gas constant. By us¬ 
ing the constitutive relations, \m = (dM/dH)r, and 
Xp = £q 1 (dP/dE)T, to define the dimensionless suscep¬ 
tibilities Xm and xp f° r the order parameters M and P 
respectively, we incorporate these variable replacements 
and accompanying constitutive relations into Eq. [T] and 
arrive at the thermodynamic relations: 

Xm(T) = p 0 T(dM(T)/dT\ H ) 2 /(c H - c M ) 


XP {T) = ^T{dP{T)/dT\ E ) 2 /{c E ^ c P ) . (2) 

Since the isochoric specific heats, Ch,m = Ch,m/V and 
Ce,p = Ce,p/V, have a very weak temperature depen¬ 
dence in the vicinity of the respective phase transition 
temperatures Tm and Tpgoj], the above expressions of 
Eq. [2] reduce to the simple proportionalities: 

X m{T) cx T M (dM(T)/dT\ H ) 2 

X p(T) cx T P {dP{T)/dT\ E ) 2 . (3) 

We emphasize that these expressions for xm{T) and 
Xp{T) are derived from thermodynamics and do not re¬ 
quire M and P to be linearly dependent on H and E 
respectively. 

B. Angle dependence of the Helmholtz free energy 

Using these same variable substitutions we write the 
total free energy as a sum of terms, 

F(P, M ) = F P (P) + F m (M) + F int (4) 

where Pp(P) and Fm{M) include quadratic terms in P 
and M respectively, i.e., 

F,.(r) + F M (M) = (fo-)p + (^)«“ (5) 

and Fi n t is an interaction term including the symmetry- 
allowed trilinear and biquadratic terms, i.e., 

Fint = ai L ■ (P x M) + a 2 (P ■ M) 2 . (6) 

where a\ and 02 are positive constants. In Eq. [5]we have 
omitted higher order quartic terms in P 4 and M 4 since, 
as will be shown below, linear response provides a good 
description of our data. 

We next minimize the angular part of the free energy 

Fi nt = a\LPM sin0 + a 2 (PM) 2 cos 2 6 . (7) 

where 6 is the angle subtended by P and M, L is rede¬ 
fined as the projection of L on P x M, and P = |P| and 
M = \M\ will henceforth be used as parametric scalar 
variables having their own temperature and field depen¬ 
dencies. The extrema of Eq. [7] at 6 = d Q are found by 
setting dFi n t/dd = 0, leading to two sets of solutions 
for 0o; namely that coming from cos0o = 0, i.e., 0$ = 
7t/2, 37t/2, or sin0o = a\L/2a2PM < 1. The free energy 
is found to be a minimum (i.e., d 2 Fi n t/dd 2 > 0) only 
for the specific cases do = 7t/2, provided a\L/2a2PM < 
1, and do = 37t/2, with no restrictions on the ratio 
a\L/2cL2PM. For the solution do = sin” 1 (aiP/2a2PM) 
with a\L/ 2 a 2 PM < 1 the free energy is a maximum. 
Accordingly, the only two possibilities giving a minimum 
in the free energy are do = 7t/2,37t/2. Substitution of 
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these angles into Eq. [7] reveals that 9q = 37 r /2 gives a 
global minimum, F int = —a\LPM , so that the total free 
energy of Eqs. [4] and [ 6 ] at 9q = 3ir/2 simplifies to 

F(P, M) = P 2 + (^ 7 ) M 2 - fi 0 aPM, ( 8 ) 

where the positive constant a = aiL/po is proportional 
to L and therefore breaks both time and inversion sym¬ 
metry. This form for the minimum free energy holds with 
no restrictions on the ratio a\L/2aiPM 1 thereby allow¬ 
ing the values of P and M to be small. 

C. Strain and a hidden antiferromagnetic order 
parameter? 

A previous treatment of BMO using Landau theory in¬ 
corporates quartic, P 4 and M 4 , and biquadratic, (PM) 2 , 
terms in the free energy [4l|. In contrast, our finding in 
the above analysis that the angle-averaged free energy 
is minimum at 9q = 3n/2 implies that the conventional 
biquadratic coupling term (P • M ) 2 in Eq. [ 6 ] is zero and 
hence can be ignored when the trilinear term L-(PxM) is 
present. Having taken account of the angular dependence 
the resulting free energy expressed in Eq. [5] is consider¬ 
ably simplified with a bilinear interaction proportional 
to PM and at the same time symmetry preserving be¬ 
cause of the absorption of L = \L\ into the constant a. 
In the following analysis we therefore consider P and M 
to be the principal order parameters with L playing an 
important but sub-dominant (“hidden” 0 ) role. 

To explicitly include strain, we recall that pristine crys¬ 
tals of BiMnC >3 are centrosymmetric IliaIl4l | and thus 
do not undergo ferroelectric transitions. In contrast our 
films exhibit a “strain-disordered” islanded morphology, 
confirmed by AFM scans 0, 0 , together with a pro¬ 
nounced ferroelectric transition (see Fig. [3]). Accordingly, 
we conclude from these observations that strain described 
by a local strain variable Sx associated with disorder, 
which is present only in films but not in crystals, is the 
source of an extrinsic mechanism for ferroelectricity in 
BiMnC >3 thin films. 

This line of reasoning is pursued by making the ansatz 
that local distortions described by Sxi couple linearly to 
site specific moments, pt and m*, to give the interac¬ 
tion term of Eq. [ 8 ] More specifically, we now write the 
interaction contribution to the free energy in terms of 
disorder-related strain terms, 

$F int = giPiSxi + g^niiSxi + k5x 2 / 2, (9) 

which when evaluated at equilibrium, i.e., Sxi = —{g'iPi+ 
92 nii)/h, are found to lower the free energy by an amount 

Fint = -kSxi 2 / 2 = -{g'lPi + favnif /2k. (10) 

For clarity, the coupling constants g\ and g '2 and the 
restoring force constant k associated with strain disor¬ 
der are primed. Expansion of the quadratic form in 


Eq. HO] and averaging over all sites gives rise to correc¬ 
tions —fa /4k and — g '2 /4fc to P 2 and M 2 respectively in 
Eq. [5] which can be ignored, and to the additional bilin¬ 
ear interaction term fafaPM / k. Using this model we can 
therefore make the identification po ol = ai L = g\ g' 2 /k in 
Eq. [8] 

The strain induced corrections to the free energy must 
preserve both inversion and time reversal symmetry. This 
is true for the g\piSxi and kbx 2 / 2 terms of our above 
ansatz but not for the magnetic coupling term, fam^Xi , 
where rrii and Sxt respectively break time reversal and 
inversion symmetry. We remedy this situation by in¬ 
cluding in our ansatz the additional assumption that g '2 
is proportional to the antiferromagnetic order parameter 
L which breaks both inversion and time reversal symme¬ 
try. Accordingly, since both a and g '2 both have the same 
symmetry properties of L, the interaction contribution, 

Pint = -poaPM = —fafaP M/2k, (11) 

includes a strain contribution which by this simple qual¬ 
itative model is symmetry preserving. 

A hidden antiferromagnetic order parameter for 
magnetoelectric BMO is postulated by Solovyev and 
Pchelkina0| using microscopic theory in which a rel¬ 
ativistic spin-orbit interaction gives rise to canted spin 
ferromagnetism. In a separate and more recent theoret¬ 
ical treatment using hybrid density functional theory, a 
highly-strained antiferromagnetic supertetragonal phase 
is predicted to have a high polarization [0 . These two 
theories of improper multiferroicity in BMO 0,0] thus 
allow for the coexistence and interplay of ferroelectric¬ 
ity and ferromagnetism in a system which is otherwise 
centrosymmetric (C2/c) and incompatible with ferroelec¬ 
tricity. 

This scenario is made more plausible with the realiza¬ 
tion that ferromagnetism in BiMnCU is due to the orbital 
ordering of the spin carrying Mn atoms that is sensitive to 
atomic spacingsjlO]. In addition the magnetization, with 
a maximum of ~ l/rs/Mn0} in our films, is significantly 
reduced from the bulk value of 3.6/Us/Mn, suggesting 
that the magnetization is inhomogeneous and likely con¬ 
tains a strain-induced antiferromagnetic component. The 
atomic spacings are modified by the high strain island 
edges 0, misaligning the magnetic moments, resulting 
in a spatially varying magnetization across the islands 
which in turn produces internal electric fields Q that are 
capable of aligning electric dipoles. There may also be 
an additional contribution to the measured magnetoelec¬ 
tric coefficient a(T) (Fig. []J from the overall reduction 
in volume of the ferroelectric regions on the application 
of a magnetic field, thus increasing the magnitude of the 
negative magnetoelectric coupling. 

Temperature dependent high-pressure neutron diffrac¬ 
tion studies have also reported a pressure induced 
monoclinic-to-monoclinic structural transition which 
leads to a transformation from a ferromagnetic to an¬ 
tiferromagnetic phaseJUJ- It is not difficult to imagine 
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that this competition and sensitivity of magnetic order 
parameters to slightly different structures, albeit at high 
pressure, is present in strain-disordered thin films at at¬ 
mospheric pressure. 

D. Bilinear interaction and enhanced 
magnetoelectric coupling 

The omission of quartic terms proportional to P 4 and 
M 4 and a biquadratic coupling term proportional to 
(PM) 2 in Eq. [8] considerably simplifies our analysis. Ap¬ 
plying the thermodynamic expressions E = ( dF/dP)x,M 
and H = (dF/dM)p,p to Eq. [HI a matrix inversion 

leads to expressions for P(E,H) and M(E,H) that are 
linear in the applied fields E and H. The result gives the 
expressions 

SP = P-P R = e 0 xp(T)E + d(T)H 

HoSM = no(M — M r ) = a.E + HoXm(T)H (12) 

where Pr and Mr are the remnant moments when the 
applied fields E and H are zero. The renormalized sus¬ 
ceptibilities xp(T) = xp/ (1 - oi 2 c- 2 xp(T)xm(T)) and 
Xm(T) = XM(T)/(l-a 2 c~ 2 xp(T)xM(T)) together with 
the magnetoelectric coeficient 

d(T) = ac- 2 Xp(T) X M(T)/(l - a 2 c- 2 XP (T)xM(T)) 

= a/(c 2 Xp\T) X ^(T)- a 2 ) (13) 

are the experimental observables. We note that P and 
a(T)H in the first line of Eq. [12] both break inversion 
symmetry since the time reversal asymmetry in a(T) and 
H appears twice and therefore cancels. A similar argu¬ 
ment applies to the a(T)E term in the second line of 
Eq. [12] which breaks time reversal symmetry with the 
preservation of inversion symmetry. 

Parenthetically, we note from Eq. [13] above that a and 
a. have the same positive sign determined by the relation 
a = aiL/iio where ai in Eq. [0] is defined as positive. If 
a\ in Eq. [S]is replaced by — a\, the free energy minimum 
occurs at Oq = 7t/ 2 rather than 9q = 37r/2 and Eq. [Hire- 
mains unchanged. The experimentally observed decrease 
in Pr with applied field H, i.e., dPp/dH = —a is thus 
consistent with the signs in Eq. [lH] 

The materials-specific coupling constant a has dimen¬ 
sions corresponding to velocity and gives rise to a linear 
magnetoelectric effect [lj . If there is no ME coupling, then 
a in Eq. [5] is equal to zero as is a(T). In this case there 
are two separate and independent phase transitions with 
the higher temperature ferroelectric transition occurring 
at T = Tp where Xp 1 = ( T/Tp — 1 )/Cp goes to zero 
and the lower temperature ferromagnetic transition oc¬ 
curring at T = T m where x~^ = (T/T M - 1 )/C M goes 
to zero. The Curie temperature prefactors, Cp and Cm, 
are taken to be dimensionless. 

In the presence of finite coupling, the essential physics 
is markedly different as captured in the behavior of 
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FIG. 5: The temperature derivatives of the magnetization 
(red) and ferroelectric polarization (black) are calculated from 
Fig. 2 and provide a measure of the magnetic and electric sus¬ 
ceptibilities (see text). The ME coupling coefficient, a(T), is 
shown in blue, and is independently calculated from the slopes 
of Pr vs. B curves at each temperature (see Fig. 4. The green 
line representing the normalized product X m(T) X p(T) (deter¬ 
mined by the independently measured temperature deriva¬ 
tives of Eqs. E]) correlates well with the measured values of 
a(T) in accordance with Eq. 1171 

the determinant (^Xp^Xm ^ ol 2 = a/a(T) associated 
with the denominator of Eq. 1131 which with substitu¬ 
tion of the expressions in the preceding paragraph for 
the temperature-dependent susceptibilities becomes 

(T - T P )(T - T m ) - a 2 CpC M T P T M /c 2 = 

(T - T* P )(T - T* m ) = aCpC M TpT M /c 2 a(T). (14) 

By factoring the left most term of the above equation into 
the product (T — Tp)(T — T^) the quadratic equation in 
T can be solved for the two roots Tp and as, 

Tp, T* m = [T p +T m ± (T P + T m ) 

(1 + 4 a 2 C P C M TpTM/c 2 (Tp - T M ) 2 ) 1/2 ]/2 (15) 

In the absence of magnetoelectric coupling, a = 0 
and the renormalized temperatures Tp and T^ are un¬ 
changed from the original separate and independent 
phase transitions respectively at T = Tp and T = Tm- 
If Tp > Tm in the uncoupled system, as is the case 
in our experiment, then as coupling increases Eq. 1151 
shows that Tp — T^ also increases with Tp > Tp and 
< Tm- A similar repulsion (increase in separation) 
occurs if Tp < Tm- It is important to recognize that 
the highest of the two temperatures, Tp, marks a sus¬ 
ceptibility transition where both interacting electric and 
magnetic dipoles exist. It does not make sense to iden¬ 
tify the lower transition at T = T^ as a phase transition 
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temperature; any remnant feature at this temperature, 
such as a susceptibility peak, is a byproduct of interac¬ 
tions which have their onset at the higher temperature 
Tp. 

Interpretation of the meaning of the renormalized tem¬ 
peratures Tp and T(/ is further clarified by solving the 
rightmost equality of Eq.[l4]for the observable a(T ), i.e., 

a(T) = ac- 2 C P C M T P T M /(T - T*)(T - T* M ) (16) 

and comparing the denominator of this result with the 
denominator of Eq. [13] The magnetoelectric coupling 
constant a(T) reaches a maximum at temperatures where 
both denominators are small, i.e., when ((?Xp Xm ~ 
a 2 ) ft 0) (Eq. [13]) or equivalently when (T — Tp)(T — 
T* m ) ft 0 (Eq. Hill) . By writing the proportionalities, 
Xp(T) oc \T — Tpl -1 and xm(T) oc |T-T^| _1 suggested 
by the discussion following Eqs. [12] Eq. [TC] becomes a 
proportionality, 

“(T 1 ) oc Xp(T)xm(T) oc |T — Tp| _1 |T — 1 _1 , (17) 

showing that the magnetoelectric coefficient a(T ) has 
the same temperature dependence as the product 
Xp(T)xm{T) of susceptibilities. 

In the final steps leading to Eq. [T7] we have replaced 
the parentheses around the temperature differences by 
absolute values, since the susceptibilities as calculated by 
Eq.[3]are proportional to well-behaved continuous deriva¬ 
tives of smooth transitions shown in Fig. [TJ We have also 
assumed that L oc \M\ so that the normalized temper¬ 
ature dependence of LM is close to that of M. Said 
in another way, the magnetization transition shown in 
Fig. [T] includes the temperature dependence of the hid¬ 
den order parameter L. 

Agreement with the overlapping temperature depen¬ 
dences predicted by Eq. [17] is shown in Fig. [5] where the 
separately measured susceptibilities, Xp(T) (solid black 
inverted triangles) and Xm(T) (solid red triangles), and 
magnetoelectric coefficient a(T ) (solid blue circles) are 
plotted with respective maxima normalized to unity. The 
susceptibility peaks, which occur at the temperatures Tp 
and T^ (Tp > T^ f ) are marked with vertical arrows and 
are seen to symmetrically straddle the peak in a(T). The 
green line representing the computed normalized product 
Xm(T)xp(T ) correlates well with the measured values of 
a(T). The data show very clearly, in accordance with the 
result of Eq.[l7]how the overlapping portions of the prox¬ 
imate susceptibility transitions amplify magnetoelectric 
coupling. If the biquadratic interaction term [oc (PM) 2 } 
alone is used in the free energy, the ME effect can also be 
linear in H but only exists below the lowest T c 00 , in 
contrast to the experimental results in Fig. [5] where a is 
finite at temperatures well above the susceptibility peak 
at T = T^p 



FIG. 6: Schematic showing evolution of polarization P (blue 
arrows) and magnetization M (red arrows) as a function of 
decreasing temperature T (vertical axis). For T > Tm, where 
Tm = 105 K is the unrenormalized bulk transition tem¬ 
perature, the rapidly fluctuating electric dipoles of the po¬ 
lar microregions of the relaxor ferroelectric dominate. As T 
decreases the dipoles increase in size and relaxation times 
(fluctuation rates) increase (decrease). At T = Tm mag¬ 
netic moments appear, which because of the trilinear cou¬ 
pling, L ■ (P x M), fluctuate at rates comparable to the re¬ 
laxor rates. As T is lowered still further towards the relaxor 
freezing temperature Tp ft 70 K, the electric and magnetic 
dipoles simultaneously maintain their coupled orientations on 
laboratory time scales at comparable transition temperatures, 
Tp ft Tp ft Tm, where the magnetoelectric coupling coeffi¬ 
cient, d(T), is maximum. As T decreases further, both P and 
M continue to increase as &(T) begins to decrease (see Figs. [T] 
an dEJ. As discussed in the text the free energy is minimized 
when the alignment between P and M is 3-7r/2; the misalign¬ 
ment shown in the figure schematically illustrates the canting 
away from 3n/2 introduced by the presence of an antiferro¬ 
magnetic component, L. 


IV. SUMMARY AND CONCLUSIONS 

In summary, our study of epitaxial thin films of 
BiMnOa has uncovered a robust ferroelectric remnant 
polarization (23/rC/cm 2 at 5 K) with a pronounced mag¬ 
netoelectric coefficient (-1.25 ns/m at 55 K), neither of 
which are present in bulk crystalline BiMnOs. The polar¬ 
ization and magnetization transitions are in close prox¬ 
imity with the temperature dependence of the product of 
susceptibilities closely following the temperature depen¬ 
dence of the magnetoelectric coefficient. This proximity 
of transitions gives a strongly enhanced linear magne¬ 
toelectric coupling. Our use of a symmetry preserving 
trilinear coupling term of the form L ■ (P x M) enables a 
thermodynamic interpretation in which, after angle av¬ 
eraging, a bilinear term proportional to PM in the free 
energy provides an excellent qualitative description of the 
data in which there is only one phase transition marked 













by simultaneous maxima in the cross susceptibility and 
magnetoelectric coefficient (Eq. [TCI). Although experi¬ 
mentally we measure the suppression of a remnant po¬ 
larization, our analysis does not require the inclusion of 
quartic or biquadratic terms in the free energy to deter¬ 
mine the remnant moments and how they are affected by 
the application of external fields. Rather, the inclusion of 
a trilinear interaction in the free energy suffices to give a 
good account of the observed large linear magnetoelectric 
coupling. 

Strain disorder clearly plays an important role in un¬ 
derstanding the multiferroicity and magnetoelectric cou¬ 
pling in BiMnOa- The ferroelectricity is of a relaxor ori¬ 
gin with relaxations due to thermally activated reorien¬ 
tation of disorder-induced polar micro regions observable 
(see Fig. 0 at temperatures greater than ss 90 K. At 
temperatures below the relaxor freezing temperature of 
s=s 70K, the remnant polarizations are stable and exhibit 
a pronounced sensitivity to the application of externally 
applied strains [HI . We find that the distortions asso¬ 
ciated with the polar micro regions are described by a 
local strain variable, which in our simple model couples 
linearly to both P and M, giving rise to the requisite 
bilinear interaction term proportional to PM in the free 
energy. These very same distortions are responsible for 
the perturbations in the orbital ordering of the spin cary- 
ing Mn atoms that are responsible for the antiferromag¬ 
netic component L. The Dzyaloshinskii-Moriya interac¬ 
tion, represented in our analysis by the symmetry con¬ 
serving interaction term L ■ (P x M), favors non-uniform 
magnetic structures that are an essential ingredient to 
the observed magnetoelectric coupling 

Our analysis leaves two open issues needing further 
research. The first of these is the absence of convinc¬ 
ing magnetic evidence (other than a suppressed satu¬ 
rated moment) for a hidden antiferromagnetic order pa¬ 
rameter as proposed for the interaction term in Eq. [6] 
and discussed in section IIII Cl The review by Belikj22| 
cites a number of papers reporting ferroelectricity in 
high resistivity BMO films Il5[ - l2ll | occurring concomi¬ 
tantly with a reduced saturated moment compared to 
the bulk value of 3.6/rs/Mn. Interestingly, a recent pub¬ 
lication on high quality pulse laser deposited BMO films 
with good stoichiometry reports high resistivity together 
with a ferromagnetic T c = 105 K and a saturated mo¬ 
ment of 3.6uB/Mn as seen in bulk, but no evidence of 
ferroelectricity J44J. This evidence, though circumstan¬ 
tial, supports our conclusion that multiferroic behavior 
of BMO films requires the presence of an antiferromag¬ 


netic component which acts to reduce the saturated mo¬ 
ment. On another experimental front, neutron scattering 
measurements could resolve an antiferromagnetic compo¬ 
nent within a ferromagnetic background, but such mea¬ 
surements are difficult because the ferroelectricity and 
large magnetoelectric effect in BMO only appears in thin 
strained films or at interfaces, not in bulk. 

A second open issue concerns the question of whether 
the proximity of the transition temperatures is acciden¬ 
tal. We have shown that in the presence of symme¬ 
try preserving trilinear coupling, there is only one phase 
transition and both Pj j and M appear at that temper¬ 
ature. Clearly a time stable magnetoelectric effect can¬ 
not be present until the temperature is low enough (i.e., 
T < Tp = 70 K to have stable dipole moments. At first 
sight there is no such restriction on the magnetic tran¬ 
sition at higher temperatures which in bulk appears at 
a Curie transition T c = 105 K. One possibility, depicted 
schematically in Fig. [Gl is that the coupling of the rapidly 
fluctuating dipole moments to magnetic moments in this 
higher temperature region (T > Tp) is sufficiently strong 
to prevent time-stable magnetization. Accordingly, in 
this very plausible scenario the transitions in Pr and M 
appear simultaneously near T = Tp, where the relaxation 
time of the ferroelectric diverges and both of the coupled 
moments are stable, as is observed experimentally. 

In conclusion, our study of a surprisingly large linear 
magnetoelectric effect in BiMnOa opens a new perspec¬ 
tive on the design of novel multiferroics that might even¬ 
tually have practical application. The combination of 
strain, magnetic ordering and relaxor ferroelectricity con¬ 
spire to bring the transition temperatures of the respec¬ 
tive ferroic order parameters into close proximity with 
each other, thereby significantly enhancing the magneto¬ 
electric coupling. In addition, our theoretical treatment 
utilizing a Dzyaloshinskii-Moriya type trilinear interac¬ 
tion avoids the use of higher order (and hence smaller) 
biquadratic interaction terms and gives a good account 
of the enhanced ME coefficient. 
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